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The strong coupling constants, 9d s dk*, 9b s bk*, 9d%dk, 9b* s bk, 9d*dk! an d 
9B*BKn where Kq, K and K\ are scalar, pseudoscalar and axial vector kaon mesons, 
respectively are calculated in the framework of three-point QCD sum rules. In par- 
ticular, the correlation functions of the considered vertices when both B(D) and 
Kq(K)(K\) mesons are off-shell are evaluated. In the case of K%, which is either 
i£i(1270) or i^i(1400), the mixing between these two states are also taken into ac- 
count. A comparison of the obtained result with the existing prediction on gD*DK 
as the only coupling constant among the considered vertices, previously calculated 
in the literature, is also made. 

PACS numbers: 11.55.Hx, 13.75.Lb, 13.25.Ft, 13.25.Hw 
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I. INTRODUCTION 

The strong coupling constants among the bottom and charmed mesons with light scalar, 
pseudoscalar and axial strange mesons are the main ingredients in analysis of their strong 
interactions. More accurate determination of these coupling constants is needed to better 
understand the strong interactions among the participated mesons, construct the strong 
potentials among them and obtain knowledge about the nature and structure of the encoun- 
tered particles. Experimentally, it is believed that in the production of the charmonium 
states like J/ip and ip' from the B c or newly discovered charmonium X, Y and Z states by 
the BABAR and BELLE collaborations, there are intermediate two body states containing 
D, D s , D* and D* s mesons (for example, the kaon can annihilate the charmonium in a nuclear 
medium to give D and D s mesons), which decay to the final J/ip and ifi' states exchanging 
one or more virtual mesons. A similar story would happen in decays of heavy bottonium. 
To exactly follow and analyze the procedure in the experiment, we need to have knowledge 
about the coupling constants among the particles involved. 

The strong coupling constants among mesons take place in low energies very far from 
the perturbative region, where the strong coupling constant between quarks and gluons 
is large and perturbation theory fails. Hence in the hadronic scale, one should consult to 
some nonperturbative methods in QCD to describe nonperturbative phenomena. Among the 
nonperturbative methods, the QCD sum rules approach [H-0] is one of the most powerful, 
applicable and attractive one as it is based on QCD Lagrangian and is free of a model 
dependent parameter. This approach has rendered many successful predictions such as 
its predictions about the vector mesons \^^. The three point correlation function has 



been widely used to calculate many parameters of hadrons (see for instance |10l-ll3j). The 
QCD sum rules for some strong coupling constants were derived by means of the three point 
functions in In the present work, we investigate various strong coupling constants among 
bottom (charmed)-bottom strange (charmed strange) mesons with scalar, pseudoscalar and 
axial vector kaons. Calculation of such coupling constants can help us in understanding the 
nature of the strong interaction among the participating particles. 

In the case of the scalar kaon, we consider the B s — B — Kq and D s — D — Kq vertices for 
both Kq(800) and _K"q(1430). Understanding the internal structure of the scalar mesons has 
been a striking issue in the last 30-40 years. Despite their investigation both theoretically and 
experimentally, most of their properties are not very clear yet. Detection and identification 
of the scalar mesons are difficult, experimentally, so the theoretical and phenomenological 
works can play a crucial role in this regard. In this work, we also calculate the coupling 
constants Qb*bk arid c/d*dk for pseudoscalar K. The next aim in the present work is to 
consider the vertices B* - B - K x and D* - D - K x for both 1^(1270) and ^(1400) axial 
states taking into account their mixture. 

Experimentally, the f^i(1270) and i^i(1400) are the mixtures of the strange members of 
two axial- vector SU(3) octets 3 Pi(K^) and 1 Pi(Kf). To avoid any confusion between the 

B meson and the sign B in the K± } we will use the instead of K^ B ^ in this article. 



The #i(1270, 1400) are related to the K"' b states via [I5i |16 



Ki(1270)) 
Ki(1400)) 



K*)sme+ 
KfjcosO- 



K\)cos6, 



(1) 
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where the mixing angle 6 takes the values in the interval 37° < 6 < 58°, —58° < 6 < —37° 
The sign ambiguity for the mixing angle is correlated with the fact that one can add 
arbitrary phase to the | Kf) and A K\). Studies on B -> ^(1270)7 and r A^(1270)z/ T 



lead to the following value for 6 [20 



# = -(34 ±13)° (2) 

In the present work, contributing the quark-quark and quark-gluon condensate diagrams 
as nonperturbative effects, we evaluate the corresponding correlation functions when both 
B(D) and Kq (K)(Ki) mesons are off-shell. Note that recently, we have investigated the 
D*DK* (892), and B*BK* (892) vertices for K* being the vector meson in the framework of 
the three-point QCD sum rules in [2l[. Moreover, the following coupling constants have been 
investigated via three-point and li ght cone QCD sum rules in the literature: D*Dir [22I. I23I. 



DDpfeg, DDJ/ip [2^D*DJ/riJffl, DW*TrWi l28j, D*D*J/ip |29|, D S D*K,D*DK [27j, 
3jJ, D D s K,D s0 DKM, DDoj MjD*D*p \m^D*Dp B sQ BK,B sl B*K [35], B s0 BK 
2, a r]7r°, a r]'7r° [3j7ao^ + ^~ B and f Q K + K~ [11,0. 

The outline of the paper is as follows. In section II, introducing responsible correlation 
functions, we obtain QCD sum rules for the strong coupling constant of the considered 
vertices. For each of the scalar, pseudoscalar and axial kaon cases, we will calculate the 
correlation function when both the B(D) and Kq(K)(K\) mesons are off-shell. In the case 
of the Ki meson, first we will calculate the QCD sum rules for the vertices B* — B — K\ and 
D* — D — K\, then using the relations in Eq. ([T]), we will acquire the QCD sum rules for the 
vertices B* — B — ^(1270) and D* — D — ^(1400). In obtaining the sum rules for physical 
quantities, we will consider both light quark-light quark and light quark-gluon condensate 
diagrams as nonperturbative contributions. Finally, in section III, we numerically analyze 
the obtained sum rules for the considered strong coupling constants. We will obtain the 
numerical values for each coupling constant when both the B(D) and Kq(K)(Ki) states are 
off-shell. Then taking the average of the two off-shell cases, we will obtain final numerical 
values for each coupling constant. In this section, we also compare our result on gn*DK with 
existing prediction in the literature. 

II. QCD SUM RULES FOR THE STRONG COUPLING CONSTANTS 

In this section, we obtain QCD sum rules for the strong coupling constants associated 
with the D S (B S ) - D{B) - K*, D*(B* S ) - D(B) - K and D*(B*) - D{B) - K x vertices. 
We start our discussion considering the sufficient correlation functions responsible for the 
corresponding strong transition involving each Kq, K and K\ mesons when both D{B) 
and K*{K){K 1) are off-shell. The following three-point correlation functions describe the 
considered strong transitions: 

• correlation functions corresponding to the D S (B S ) — D(B) — Kq vertex: 

- for D{B) off-shell: 

U D(B) =l 2 J d 4 x d 4 y e i P '. x e iq-v( Q \ r (jf5( x ) V D( - B \y) r) D ^(0)) |0>, (3) 

- for Kq off-shell: 

n x o = « 2 y d* x d 4 y e ip '- x e iq - y (0\T(ri D{B) (x) rj K °(y) r/ D «™(0)) |0), (4) 
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• correlation functions corresponding to the D*(B*) — D(B) — K vertex: 

- for D(B) off-shell: 

u d(b) =i 2 j d 4 x d A y eip '- x e iQ y(Q\ r (rj K ( x ) r} D ^ B \y) ^ :(B;)t (0)) |0>, (5) 

- for K off-shell: 

njf = i 2 J d A x d A y e ip '- x e iq - y (0\T(r] D(B) (x) r] K (y) rffX B ™(oj) |0>, (6) 

• correlation functions corresponding to the D*(B*) — D(B) — K 1 vertex: 

- for D(B) off-shell: 

= i 2 J d A x d A y e ip '- x e^iOlT^u 1 ^) V D(B) (y) < s * (Bs * }t (0)) |0), 

(7) 

<j? = ?j d 4 V ^ e iq V (0\T(v^) V D(B) (y) rj^(0)) |0>, 

(8) 

- for K x off-shell: 

n^ 1 =t 2 J d A x d A y e ip '- x e iq - y (0\T(r] D(B) (x) rtf^y) <* (B * }t (0)) |0), 

(9) 

K™> = * 2 / * x d 'y e¥ ' X e iq - y (0\T(v D{B) (x) rfc{y) *£ :(B;)t (0)) |0), 

(10) 

where T is the time ordering product, q is the momentum of the off-shell state, p' is the 
momentum of the final on-shell state. We will set the momentum of the initial state as 
p = p' + q. In the vertex containing K 1 meson, we have two correlation functions for both 
off-shell cases since this meson couples into two interpolating currents v^ 1 and r/^,. We 
will define these couplings in terms of G-parity conserving and G-parity violating decay 
constants later. The interpolating currents, which produce the considered mesons from the 
vacuum with the same quantum numbers as the interpolating currents can be written in 
terms of the quark field operators as following form: 



rj K ° (x) 


= s(x)Uu(x), 


r] K (x) 


= s(x)~f 5 u(x), 


rj^(x) 


= s(x)-f u -f 5 u(x), 


rfc{x) 


= s(x)a vv f^ b u(x) 


V D ^ B \x) 


= u(x)~f 5 c(b)(x), 


T] D ^\x) 


= s(x)-f 5 c(b)(x), 


rj^(x) 


= s(x)^c(b)(x), 



where U stands for unit matrix and u, s, c and b are the up, strange, charm and bottom 
quark fields, respectively. 
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According to general philosophy of the QCD sum rules, we calculate the aforementioned 
correlation functions in two different representations. In physical or phenomenological rep- 
resentation, we calculate them in terms of hadronic parameters. In QCD or theoretical 
representation, we evaluate them in terms of QCD degrees of freedom like quark masses, 
quark condensates, etc. with the help of the operator product expansion (OPE), where 
the perturbative and nonperturbative contributions are separated. The QCD sum rules for 
strong coupling constants are obtained equating these two different representations through 
dispersion relation. To suppress contributions of the higher states and continuum, we will 
apply double Borel transformation with respect to the momentum squared of the initial and 
final on-shell states to both sides of the obtained sum rules. 

First, let us focus on the calculation of the physical sides of the aforesaid correlation 
functions for example when D(B) meson is off-shell. Saturating the correlation functions 
with the complete sets of three participating particles and isolating the ground states and 
after some straightforward calculations, we obtain: 

• physical representation corresponding to the D S (B S ) — D(B) — Kq vertex: 

(q 2 - m 2 D(B) )(p 2 - m 2 Ds(Bs) )(p' 2 - m 2 K *) 
+ ..., (12) 

• physical representation corresponding to the D*(B*) — D(B) — K vertex: 
j-j--D(b) = (Q\rnK{p^ 

+ ..., (13) 

• physical representation corresponding to the D*(B*) — D(B) — vertex: 

uD{B ) = (oi^ (B) i^(i?)(g))(^:(i?:)(p,6)i^^)|o) 

^ i<f - m l(B))(P 2 ~ m D*(B- a )) 



X 



+ 



rj^K-(p\e'))(K-(p' : e')D(B)(q)\D*(B*)(p : e)) 

+ ..., (14) 



(p' 2 -m 2 K a) 

(0|^|^( P ', e '))(^(p , ^ , ) J D(5)(?)l^(5:)(p,e)) 



(p' 2 - m 2 Kb ) 



n 



D (B) _ (0|^ fl )|D(B)( g )>(D:(B:)(p,e)|^:W)|0> 



(q 2 - m 2 D{B) )(p 2 - m 2 D . AB:) ) 



x 



+ 



(p' 2 -m 2 K a) 

(0\vM(p',^)(KKp',elD(B)(q)\D:(B:)(p,e)) 



(p' 2 - m 2 R ,) 



+ 



(15) 
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where .... represents the contributions of the higher states and continuum, and e and e' are 
the polarization vectors associated with the D* and mesons, respectively. From the 

above equations it is clear that to proceed we need to define the following matrix elements 
in terms of decay constants as well as strong coupling constants: 



(0\r] K °\K* Q (p')) = m K5 f K *., 
(0\V K \K( P ')) - mlfK 



(0\r, D ^\D s (B s )(q)) = fD ° m 

m c{ b) + m s 

(£>:(s;)(p,e)|^W)|0> = m Dt{Bt) f Dt{Bt) 
(K* (p')D(B)(q)\D s (B s )(p)) = g DsDK * {BaBK *)P ■ p' 



(K(p')D(B)(q)\D* s (B* s )(p, e)> = g DtDK{ B t BK){p - q) - e, 
(Kf\p',e')D(B)(q)\D:(B:)(p,e)) = ^ J ,^« (fl ; BX? w ) {m^ ( fl :) (c*.c) + (e* .p)(e.p% 

(16) 



and 



(0\7]^\K^p', e')) = m K af K ae' u , ) 

(0\v%\KW,e>)) = /^(eX - <X) J C ° nSerVing definitions ' 

(0\v«M*(p', e')> = f K ? 4> K *^ v p' v , - e>'J, . r , . , r , r r 

„ ii K b > G-parity violating definitions, 

{^\K\{p'^')) = m Kb J Kbi ^ K ' ( 



l e u 



(17) 



where f K *, f K , f D ( B ), /d s (b s ) and /d»(b*) are leptonic decay constants of the K%, K, D(B), 
D S (B S ) and D*(B*) mesons, respectively. The fx* and f K b± are decay constants encountered 
to the calculations from both definitions of the G-parity conserving and violating matrix 
elements for the axial and K\ states (for more details see 15, H(J 0, SH). The % ,Kl 

II K b 

and %' 1 are zeroth order Gegenbauer moments. In the above equations, the gD s DK*(B s BK*), 
9d*dk(b*bk) and 9 D * DK *ib), B * BK a(b}} are strong coupling constants, which we are going to 
obtain QCD sum rules for them in this section. 

Using Eqs. (|T2|) - (fl6|) . and summing over the polarization vectors using the 

j J* _ , PvP'e 



e u e e — ~9ue + 



m K a(b) 
^1 

* __ I PvPl* /I Q\ 

m Di{B* s ) 

the final physical representations of the correlation functions for each vertices in the case of 
D(B) off-shell is obtained as: 
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D a (B s ) — D(B) — K* vertex: 



Y\D(B) 



MB) 



J K* K* mc ^ b j+ mu m c ( 6 )+m s 



+ m 



A'* 



J(B) 



D*(B*) - D(B) - K vertex: 



(19) 



n I>(B) _ D(B) ( 2) -_ 

LL H — yD*DK(B*BK) J ^2 _ 



/• fl<m, 2 K fp(B)m 2 D(B) 

jD*(B*)m D *( B *) 



m l(B))(P' ~ ™ 2 K )(P 2 ~ m %{B* a ) 



1 + - 2 P ; 



- D(B) - K x vertex: 

u d(b) _ ( d(b) , 2) rn K af K a D(B) 

X V ~~ \ yD* s DK^{B* s BK1)\H ) ^,2 _ m ^ ^ J DK\ (B* BK\ 

l /d;(b;)"1£»:(b;) 



II a 6 



x 



fp{B)m D{B) 
rn c(b) +m u 



(<1 2 ~ m l(B))(P 2 - m D t {B t )J 



m 



(20) 



d*(b*)9v + other structures > + (21) 



i V^' ~~ ^DJDKfOBJBtff)^ )j^—^2^ +9 D*DK b 1 (B*BK b 1 ) 



M B ) 



D(B) 



M 2 



f 



K b 



,/2 



X 



fp(B)m 2 D{B) 



f D *{B*)m D *( B *) 



(q 2 - m 2 D{B) )(j> 2 - m 



M(B*)J 



m D*(B*)(9^P'u' - 9ixv'P' v ) + other structuresj + 

(22) 



where to calculate the coupling constants, we will choose the structures p^ and g^v^iivV'v 1 ~~ 
9/j.v'Pl) from both sides of the correlation functions corresponding to the vertices containing 
the K and K x with current v^ x (K 1 with current rj^,), respectively. From the similar way, 
one can easily obtain the physical representations of the correlation functions associated 
with the i£J (#)(#!) off-shell. 

Now, we concentrate to calculate the QCD or theoretical sides of the considered correla- 
tion functions. The QCD representations of the correlation functions are obtained in deep 
Euclidean region, where p 2 — > — oo and p' 2 — > — oo via OPE. For this aim, each correlation 
function in QCD side is written in terms of the perturbative and non-perturbative parts as: 



IT = IT + IT 

per 1 nonper 



(23) 

where % stands for D{B) or Kq(K)(Ki) and the perturbative parts are defined in terms of 
double dispersion integral as following: 



n 



per 



ds 



(s — p 2 )(s' — p' 2 ) 



+ subtraction terms, 



(24) 
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FIG. 1. Diagrams considered in the calculations. The first and third line diagrams refer to the 
B(D) off-shell and the second and fourth line diagrams show the case when Kq (K ) (K\ ) is off-shell. 

where p l (s,s',q 2 ) are called spectral densities. In order to obtain the spectral density, we 
need to calculate the bare loop diagrams (a) and (d) in Fig. for D(B) and Kq(K)(Ki) 
off-shell, respectively. We calculate these diagrams in terms of the usual Feynman integrals 
by the help of the Cutkosky rules, where the quark propagators are replaced by Dirac delta 
function, i.e., 2 ^ m 2 — > ( — 27u)5(g 2 — m 2 ). As a result, the spectral densities are obtained as 
follows: 



D S (B S ) — D{B) — K* vertex: 
- D(B) off-shell: 



p DiB) (s,s',q 2 ) 



|m s (m u (m s + m u ) - g 2 ) - sm u 



2 A 1 /2( S;S ' ; g2) 

m c (b) ((m s + m u ) 2 - s' - m c{b) (m s + m„)) } , (25) 



K* off-shell: 



p K ° (s, s', q 2 ) = 2 A i/2^y >g 2) { m * s ' + m < b ) (( m * + m «) 2 _ 



m c(6)( m « + m «) +m u (s- m s (m s + m u )) \ , (26) 



D*(B*) - D(B) - K vertex: 
- D(B) off-shell: 



p D(B \s,s',q 2 ) 



W 2 {s,s',q 2 ) 



K - m s )(q 2 - s)(m c{b) m 2 s + m u (s - m 2 s - q 2 )) 
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- s'( - m^m u + 2ml {b) (m u - m s ) - 2m 2 s q 2 + m 2 c(b) (2m s m u + q 2 - s) 
+ <l 2 (s - Q 2 ) + m s m u (s + q 2 ) + m c{b) {m s - m u )(m 2 s + q 2 + s) 



I 1 1 2 

— s (m 



1(b) - m c(b)m s + m c(b) m u + q 2 ) 



(27) 



K off-shell: 

p K (s,s',q 2 ) = 



, L ( m c(b) - m u )(q 2 - s)(m 2 c(b) (m c{b ) - m u ) + m u (-m s m u - q' 
a [s, s , q ) 

+ s)) + (mj t) (m ( - m u ) + 2m 3 s m u + m 2 c(b) (-m s m u - 2q 2 ) + m 2 (g 2 - s) 
+ <? 2 (s - <? 2 ) - m s m u (q 2 + s) + m c ( b )(-2m 3 s + 2m 2 s m u + m u (g 2 + s) 



+ m s (g 2 + s)))s' + {-m c{b) m s + m 2 + m s m u + g 2 )s'" 
D*(B*) — D(B) — K x vertex: 
- D(B) off-shell: 

Pi (B \s, s', q 2 ) = -2NJ (s, s', q 2 ) 2m 3 s + 2m 2 m u - 2m c ^m s (m s + m u 
+ m s (s + s' - g 2 ) + AA{m c{ p) - m^j 
+ B(2m u s + m c(6 )(g 2 — s — s') + m s (-q 2 + 3s + s')) 
+ C( - 2m c(6) s' + m u (-g 2 + s + s') + m s (-q 2 + s + 3s')) 



(28) 



(29) 



p? ( B) (s, s', g 2 ) = 4iV c / (s, s', g 2 ) [24 + (m c(6) m s - m 2 ) - (B + 2G)s + #(g 2 - s - s') 



+ C(m c ( b) (m s + m u ) - m s (m s + m u ) + g 2 - s - s' 



Ifi off-shell: 

pf 1 (s, s', g 2 ) = 2NJ (s, s',q 2 ) 2m c(b) (m 2 c{b) + m s m u - m c{b) (m, + m u ) + s') 
- AD{m c{b ) - m u ) + E(2m c ( b ) - m s - m u )(q 2 - s - s') 



(30) 



+ 2Fs / (-2m c(6) + m s + m u ) 



pf ^s, s', g 2 ) = AN c I (s, s', g 2 ) [ - 2D + m 2 (6) - m c(6) m s - £s 



- F(m c ( b )(m c (p) -m s - m u ) + m s m u 
where pi and p 2 correspond to the currents rj^ 1 and r/^),, respectively and 



(31) 



(32) 



-m 2 g 2 (-g 2 + s + s')}, 
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l{m 2 (V - s + s') - s'(2m 2 (b) - g 2 - s + s')}, 



10 



C = ^{™c(&)( ~ q 2 + s + s') + s(q 2 - s + s') - m*( - q 2 - s + a')}, 

n lf4 2 i 4 / i 2 / 2 2/ 2 , , /\ . 2/2 / 2, /\ 

2A^ mc{b)g +m « s+?ss_m c(fe)? (-9 + s + s ) + m s {m c{b) {-q +s-s) 
+ s'(-q 2 -s + s'))}, 



E = ^{s'(2m 2 - q 2 - s + s ')- m 2 c(b) (q 2 - s + s% 

F = T77V^ ( ~ m ' + + S ^~ q2 + s + s ')- 2s ( m l(b) + sf)h 
Ays, s , q ) 

G = -^{3m\ b) s'(q 2 -s-s') + m 4 s (2q 4 - (s - s') 2 - q 2 {s + s')) -ss'(- 2q 4 + (s - s') 2 
+ Q 2 (s + s'j) + m 2 s (q 6 - q\s + s') + (s - s') 2 (s + s') - q 2 (s 2 - 6ss' + s' 2 )) - 2m\ b) 
( s'(q 4 - 2s 2 + q 2 (s - 2s') + ss 1 + s' 2 ) + m 2 s (q 4 + s 2 + ss' - 2s' 2 + q 2 (-2s + s')))}, 

H = ^rf^ ~ 2 <? 2 (* ~ 2s ') + ( s ~ s ') 2 ) + <{q A + 2g 2 (2. - s') + (s- s') 2 ) + mt (b) 

( q 4 + s 2 + Ass' + s' 2 - 2q 2 (s + s')) - 2m 2 s s( - 2g 4 + (s - s') 2 + q 2 (s + s')) - 2m\ b) 
( m 2 (q 4 - 2s 2 + q 2 (s - s') + ss' + s' 2 ) + s(q 4 + s 2 + ss' - 2s' 2 + q 2 (-2s + s')))}, (33) 

and 

/J* 1 



Io(s,s,q 



4A 1 /2( S)S ' ; g2) 

A = q 4 + (s - s') 2 - 2q 2 (s + s') 
A(a, b, c) = a 2 + b 2 + c 2 - 2ac - 26c - 2ab. (34) 

In the spectral densities, iV c = 3 is the color number and we have kept terms linear in m u . 

Now, we proceed to calculate the nonperturbative contributions in QCD side. We consider 
the quark-quark and quark-gluon condensate diagrams presented as (b), (c), (e), (f), (g), 
(h), (i), (j), (k), (1), (m) and (n) in Fig. (OQ). Contributions of the diagrams (c), (e), (f), 
(g), (i), (k), (1), (m) and (n) in Fig. ([!]) are zero since applying double Borel transformation 
with respect to the both variables p 2 and p' will kill their contributions because of only one 
variable appearing in the denominator in these cases. Hence, we consider contributions of 
only diagrams (b), (h) and (j) in Fig. fll}). As a result, we obtain: 

• D S (B S ) — D(B) — K* vertex: 

- D(B) off-shell: 

n D(B) = (gg) ( 2 rn c{b) m u -m 2 c{b) + q 2 _ 1 _ 1 m 2 (Am c{b) m u - m 2 c{b) + q 2 ) 

nonper ~ ^ I rr ' r r ' 



ml ml ml(m 2 c{b) - Am c{h) m u - q 2 ) m , 



,2 ™2 



K* off-shell: 



Sniper = 0, (36) 



D* S (B*) - D(B) - K vertex: 
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- D(B) off-shell: 



nonper 



- if off-shell: 



n 



K 

nonper 



o, 



2rr r 



D*(B* S ) — D(B) — Ki vertex: 
- D(B) off-shell: 

tjD(B) 
11 nonper 1 

q 2 + 2m u m c{b) - m 2 c{b) 



2 i ,™ ... »^.2 



+ 



2rr' 



2 


+ 1 


rr' 




1 


l-i 


27 





r /2 



1 1 
m e(b) - T 



(37) 



(38) 



(39) 



11 nonper2 



-(ss) 



m c(6) m5 m c(b) 



rr' 



24 rr' 2 



ifi off-shell: 



n ^i 

nonperl 



o, 



n Ki = n 

nonper2 "i 



(40) 

(41) 
(42) 



where r = p 2 — m 2 ^ and r' = p' 2 — m 2 . The H n0 nperi and n nonper2 correspond to also the 
currents r/^ 1 and 77^,, respectively. In this step, we equate the physical side in the case 
of Kq and the coefficients of the selected structures in physical sides of K and K\ to the 
corresponding QCD sides and obtain QCD sum rules for the considered strong coupling 
constants. To suppress the contributions of the higher states and continuum, we also apply 
the double Borel transformation with respect to the variables p 2 (p 2 — > M 2 ) and p' (p' — > 
M' ). Finally, we get the following sum rules for the considered coupling constants: 

• D S (B S ) - D(B) - K* vertex: 



D(B) off-shell: 

MB) 



9D s DKl(B s BKl)yl 
1 



2(q 2 - m 2 D{B) )(rn c (p) + m u )(m c[b) + m s ) 
m ls(B a ) m2 D(B) m K* fD s (B s )fD(B)fK*(m 2 Ds(Bs) + m 2 K * - q 2 ) 



Ds(Bs) 



e M 



2 



x e m 



+ BU D ^ 

1 nonper 



so 



ds 



47T 2 J \m c , b) +m s ) 2 J(m s +m u ) 2 



[ S '° , ds'p D ^ B \s,s',q 2 )9[l - U D(B \s,s')) 2 ]e^e^ 



(43) 
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- K* off-shell: 



1 



-2(g 2 - m |,)K (6) + mM )K (t)+ms 



m k(B a ) m l(B) m K*jD s (B s )fD(B)fK*(m 2 Ds{Bs) + m 2 D{B) - q 1 



Da(Bs) 

-e 



x e m 



2 

Ps(-Bs) 
72 



. ■-, , ds . 

47T J (m c{b) +m s )' 2 J (m c{b) +m u ) 



ds'p K ~o{s,s',q 2 )6[l - (f K o(s,s')) 2 ]e^e^ 



• D*(B* S ) - D(B) - K vertex: 
- D(B) off-shell: 

d(b) , 2 , _ (<f - ml {B) )(m c{b) + m u )(m s + m u 

9D* a DK(B*BK)\1 ) 



/dj {b* ) Jd(b) fKm D * (B*)m 2 K m 2 D( B) ( 1 + 



9 9 

n D*(B*) 



m2 D*s[B* a ) m K 



x 



K off-shell: 



9d*dk(b*bk) id 2 ) 



l — r ds [ S '° ds'p D ^ B \s,s f ,q 2 )9[l - {f D{ - B \s,s')f}e^e^ 

(45) 



(q 2 - m 2 K )(m ci b) + m u )(m s + m u ) 



2 2 

m g|(B|) m g(B) 

-e e m' 2 



/d; (-B| )fD(B)fKm Dt{Bi) m 2 K m 2 D{B) {\ + 



D(B)~ 



-D*(B*) 



X 



J (m nlh \+m s ) 2 J(m r 



4 7T J (m c{b) +m s ) 2 J(m c{b) +m u ) 2 



ds'p K {s,s',q 2 ) 9[1 - {f K {s,s')) }e^e^ 



(46) 



(44) 



£>;(£*) -D(B)- K x vertex: 
- off-shell: 



9d*^)K^ (B* BKl) ( ? 2 ) e _ ^ 7!jL + 9^DKl{B*BKl)^) m KlfK^ «0 

l £)(B)y 



(9 2 -^D(m) 



-e m 



O s *(B*) 
2 



/d; (b;)/d(b) mc( J2L m i>; (b; ) 



X 



so 



(is 



4 7T J(m c{i) )+m, s ) 2 J (m s +m u ) 2 



ds'p^ B \s,s', q 2 )e[l-(f D ( B \s,s>)) 2 ] 



x e^e^ 77 + BU n ^ erl 



(47) 
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D{B) / 2\ f -J-i-^i -pf i U \P) I 2\ f 

(? 2 -™D(B)) 



V 2 



£>(B) 



A/' 



-e 



/d; (B; ) /d(b) mc(6)+mu ^Dj (B| ) 

1 



X 



, ... ° d^( S , S ',, 2 )^[l-(/^)( S , S ')) 2 ] 

*± 7T J(m c (,,)+m s )^ J (m s +m u y 



x eFes? + .BIT 



D(B) 
nonper2 



(48) 



- fTi off-shell: 



^5b^(B ? B^)(? 2 ) (g2 ) + ^^(BJB^)(? 2 ) (g 2 _ ^ 

1 



-, K 1 



m K b 1 fK b 1 ±a 



K\> 



/d; (b; ) /d(b) ^f^T m D* (b* 



2 2 
m D*(B|) m g(B) 

-e m 2 e m' 2 



x 



— / rfs / rfs'pf 1 (*, s', " (f Kl (s, s')f] 



x e^e^ + BU^ perl 



K 



(49) 



(a 2 ) /y? Q ° I a K * (a 2 ) ^ 

jD*DK^(B*BK^) \M ) ^ 2 _ m 2. o ) ^ yDiDK\{BiBK\)Vl > ( g 2 _ m 2_ ) 



/ D* (B* )fD{B) m c( ^+ mu m Di (B* ) 
1 



X 



4 7T J(m c ( b) +rn s ) 2 J(m c ( b) +rn u ) 2 



f s 
J (m„, 



ds>p^{s,s',q 2 )0[l-{f K \s,s')f\ 



x e^eM' 2 +BU^ per2 



(50) 



where the BU nonper represents the double Borel transformation of the non-perturbative part 
in each case, s and s' Q are the continuum thresholds and the functions f l (s,s') inside the 
step functions are determined requiring that the arguments of the three S functions coming 
from the Cutkosky rule vanish simultaneously. As a result, we find: 

• D(B) off-shell: 



f D{B) (s,s') = 



K^K)^) off-shell: 



2 s (m 2 -ml + s') + (m^ - m 2 - s)(-q 2 + s + s') 



A1/2 ( m c(b)' m s>s)^ 1/2 (s>s / ,g 2 ) 



(51) 



A 



s,s') 



2 s {-m\ b) +m 2 u - s') + (m 2 c(b) - m 2 + s)(-q 2 + s + s') 



(b) 



W 2 (m 2 c{b) ,m 2 ,s)\i/ 2 (s,s>,q 2 ) 
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(52) 

Here, we should stress that the physical regions are imposed by the limits on the integrals 
and step functions in the integrands in the sum rules expressions. In order to subtract the 
contributions of the higher states and continuum, the quark-hadron duality assumption in 
the following form is used: 

p hi 9 herstates^ ^ = p OPE^ _ _ (53) 

The double Borel transformation used in calculations is also defined in the following way: 

B \ \ v (_l)m+n_}: l _ -mf/M 2 -m^/M' 2 1 (M) 

(p 2 - m\) m {p> 2 - mlY K ' T{m)T{n) (M 2 ) m - 1 (M' 2 ) n - 1 ' 1 ' 

At the end of this section, we would like to mention that using Eqs. ([T]) and ( I16p . the 
couplings to i^i(1270) and f^i(1400) are obtained in terms of the couplings to the as: 

9D*DK 1 {1270)(B^BK 1 (1270)) = 9D*DK?(B;BK<1) SU16 1 + gD*DKl(B* s BKl) cos $ 
5 , D|D_fsTi(1400)(B|BA"i(1400)) = 9D*DK?{B*BKf) COS 9 — go'DR^B^R*) s ^ n ^ (55) 

III. NUMERICAL ANALYSIS 

In the present section, we numerically analyze the expressions of QCD sum rules obtained 
for the considered strong coupling constants. Some input parameters used in the calculations 
are: m K = (493.677 ± 0.016) MeV, m K *(800) = (672 ± 40) MeV, m K *(U30) = (1425 ± 
50) MeV, ra Xl (1270) = (1272 ± 7) MeV, m^(1400) = (1403 ± 7) MeV, m D = (1.8648 ± 
0.00014) GeV, m B = (5.2792 ± 0.0003) GeV, m Ds = (1.96847 ± 0.00033) MeV, m B = 
(5.3663 ±0.0006) MeV, m D * = (2.1123 ± 0.0005)Gey, m B * = (5.4154 ± 0.0014) GeV 



m c = 1.3 GeV, m b = 4.7 GeV, m s = 0.14 Ge^|43j, f K = 160 MeV[U], /^ (8 oo)(l GeV) = 
(340 ±20) MeV, f K xqi30)(l GeV) = (445 ±50) MeV gj, f D . = (272 ±16^1) MeV, f B * = 
(229±20i\ 6 ) MeV f B = 190±13 MeV 0, f D = (202±41±17) MeV Q, f Ds = (286± 



44 ±41) MeV [49^ f Bs = 196 MeV 50], (ss) = 0.8(uu) = -0.8(0.24 ± 0.01) 3 GeV 3 |43j 
m 2 = (0.8 ± 0.2) GeV 2 0, m K a = 1.31 GeV, m K b = 1.34 GeV, f K «{l GeV) = 0.25 GeV, 

f Kb ±(J GeV) = 0.19 GeV, a5'^(l GeV) = -0.19 ± 0.07 and 4' K ' {1 GeV) = 0.27±g;?? 



|20|,|4C 

The sum rule for strong coupling constants also contain four auxiliary parameters, namely 
the continuum thresholds So and s' related to the initial and final channels, respectively as 
well as Borel mass parameters M 2 and M' 2 . These quantities are mathematical objects, so 
according to the general criteria and standard procedure in QCD sum rules, our physical 
results should be insensitive to them. Therefore, we shall look for working regions of these 
quantities at which the dependence of coupling constants on these auxiliary parameters are 
weak. The working regions for the Borel mass parameters M 2 and M' are determined 
demanding that both the contributions of the higher states and continuum are sufficiently 
suppressed and the contributions coming from the higher dimensional operators are small. 
Our calculations lead to the following working regions common for all cases: 

• D^DK^{K){K X ) vertex: 
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- D off-shell: 8GeV 2 < M 2 < 25GeV 2 and 5GeV 2 < M' 2 < 15GeV 2 , 

- Kq(K)(Ki) off-shell: QGeV 2 < M 2 < 15GeV 2 and AGeV 2 < M' 2 < 12GeV 2 , 

• B^BKq(K)(Ki) vertex: 

- B off-shell: UGeV 2 < M 2 < 30GeV 2 and 5GeV 2 < M' 2 < 20GeV 2 , 

- K*(K)K 1 off-shell: QGeV 2 < M 2 < 20GeV 2 and 5GeV 2 < M' 2 < IhGeV 2 . 

The continuum thresholds, So and s' are not completely arbitrary but they are correlated 
to the energy of the first excited states with the same quantum numbers as the considered 
interpolating currents. Our numerical calculations show that in the regions (m« + 0.3) 2 < 
so < (rrii + 0.7) 2 and (m/ + 0.3) 2 < s' < (rrif + 0.7) 2 , respectively for the continuum 
thresholds s and s' , our results have weak dependence on these parameters. Here, m, is 
the mass of initial particle and the rrif stands for the mass of the final on-shell state. For 
instance consider the g£> 3 DK* coupling constant at which the D meson is off-shell. This 
coupling constant describe the strong transition D s — > DKq, and for this case rrii = mr> s 
and ??!/ = m,K*- 

As an example, we present the dependence of strong coupling constant ^^^'(soo) 011 
Borel mass parameters at Q 2 = 1 GeV 2 , where Q 2 = —q 2 in Fig. ([2]). This figure demon- 
strates a good stability of the results with respect to the variations of Borel mass parameters 
in their working regions. 



9 12 15 18 21 24 

Total Contribution 

Perturbative Contribution 

Non-Perturbative Contribution 



1.0- 



0.0- 



1.0 



0.0 



9 12 15 18 21 24 



1.0- 



0.0- 



12 

- Total Contribution 

■ Perturbative Contribution 

■ Non-Perturbative Contribution ' 



1.0 



0.0 



FIG. 2. Left: 9d s ' dk * {800) {Q 2 = 1 GeV 2 ) as a function of the Borel mass M 2 with M' 2 = 10 GeV 2 . 
Right: 2^b X * (80 o)(<2 2 = 1 GeV 2 ) as a function of the Borel mass M' 2 with M 2 = 17 GeV 2 . The 
continuum thresholds so = 6.09 GeV 2 and s' = 1.37 GeV 2 have been used. 

Now, we proceed to find the Q 2 behavior of the considered strong coupling constants 
using the working regions for auxiliary parameters. First, we consider the scalar kaon case 
for both ^(800) and Kq(1430). The strong coupling constant in this case obeys from the 
following Boltzmann function: 

g(Q 2 ) = A, + [GeV- 1 ]. (56) 

l + exp^ 
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The values of the parameters Ax, A%, Xq and Ax for the considered coupling constant form 
factors are given in Table [B The coupling constants are defined as the values of the form 





-^1 






Ax 


9d 3 DK*[800)W > 


Q /IRS 


9 7/11 
— Z.I 41 


S DR7 
O.UD / 


4.y»o 


(iC *(800)) (n2) 
9d s DK* (800) W > 


n no/i 

— U.UZ4 


n 779 


K 79 "i 
0. / Zo 


1 9^7 
L.ZO 1 


yD s DK*(1430)^ > 


4.712 


-3.818 


24.863 


10.985 


(«S(1430)) (n2] 
9 D S DK* (1430) W 1 


-0.022 


0.772 


4.729 


1.637 


a [B) (O 2 ) 
9 B S BK* (800) V* > 


4.151 


-1.932 


13.842 


12.149 


(K*(800)) (n2) 
9 B S BK* (800)^ > 


-0.017 


0.547 


5.431 


1.121 


a (B) (O 2 ) 
VBsBK* (14Z0) > 


2.055 


-0.207 


11.239 


5.084 


(1^(1430)) (r)2] 
UBsBK* (1430) W 1 


-0.004 


0.255 


4.819 


1.146 



TABLE I. Parameters appearing in the fit function of the coupling constants for D s DKq(800), 
D S DK$(U30), B S BK%{800) and ^£1^(1430) vertices. Ax and A 2 are in GeV~ l units, while x 
and Ax are in the units of GeV 2 . 

factors at Q 2 = — wim eson , where m meson is the mass of off-shell meson. The results of the 
coupling constants obtained using Q 2 = — m 2 meson are given in Tables [TT] and IHII The final 
result for each coupling constant is obtained taking the average of the coupling constants 
obtained from two different off-shell cases. The errors in the numerical results are due to the 
uncertainties in determination of the working regions for the auxiliary parameters as well as 
the errors in the input parameters. 





Q 2 = -m 2 D 


Q 2 — m 2 K ^ S0Q ^ 


Average 


9 D a DK* (800) 


0.97 ±0.02 


0.74 ± 0.05 


0.85 ±0.08 




Q 2 = -m 2 D 


Q 2 = — m 2 K *( U3Q j 


Average 


9 'D S DK* (1430) 


1.16 ±0.12 


0.49 ± 0.07 


0.83 ± 0.09 



TABLE II. Value of the 9d s dk^(800,1430) coupling constant in GeV 1 unit. 

In the case of pseudoscalar kaon and D off-shell, the strong coupling constant is well 
described by the following monopolar fit parametrization: 

m - 8 ' 76 {GeV2) (57) 

2 _ ™2 

meson 



The value of coupling constant obtained at Q = —m meson is presented in Table fjIV 
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— —m B 


ra — ,™2 

~ ,n K*(800) 


Average 


9 B S BK* (800) 


2.28 ±0.18 


0.53 ±0.09 


1.41 ±0.21 




Q 2 = -ml 


Q 2 = — »7ljf*( 1 43o) 


Average 


9B S BK*(1AZ0) 


1.85 ±0.53 


0.25 ±0.04 


1.05 ±0.32 



TABLE III. Value of the 5 , _b s b_k'*(800,1430) coupling constant in GeV 1 unit. 





Q 2 = -ml 


Q 2 = -ml 


Average 


9d*dk (Present work) 


2.79 ±0.24 


2.99 ±0.26 


2.89 ±0.25 


9D t DK ([30]) 


2.72 


2.87 


2.84 ±0.31 



TABLE IV. Value of the go* dk coupling constant. 



The result for strong coupling constant of pseudoscalar case and an off-shell K meson 
can be well fitted by the exponential parametrization 



9d*dk 



(Q 2 ) 



3.55 e 7 - 25 (Gev^) _ o.88, 



(58) 



where using Q 2 = —ml, we obtain the result as also presented in the Table (lIVp . We also 
depict the final result for this case taking the average of two above obtained values. This 
Table also shows the predictions of [30J on g D * DK as the only existing previously calculated 



coupling constant among the considered vertices. Comparing our results with that of [30 
we see a good consistency between two predictions. 

Similarly, for B*BK vertex, our result for the pseudoscalar kaon and B off-shell is better 
extrapolated by the exponential fit parametrization 



9b*bk(Q 2 ) = 0-66 e + 0.23, 

and in the K off-shell case by the parametrization 

Q 2 

4.39 e 4-02 (Gev z ) _ i.03. 



9b*bk 



(Q 2 ) 



(59) 



(60) 



Using the same procedure as above, we find the values depicted in the Table (jvjl . 

In the case of axial vector kaon, the strong coupling constant obey also the same Boltz- 
mann function as the scalar case. The values of the parameters A\, A 2 , Xq and Ax for 
coupling constants in this case are given in Table IVI1 The same procedure as in the scalar 
and pseudoscalar cases leads to the numerical results for the corresponding coupling con- 
stants as presented in the Tables flVIII and IVIIip . 

In summary, the strong coupling constants, gB a BK*, 9d s dk*, 9b*bk, 9d*dk > 9b%bk x and 
9 d* DK\ i have been calculated in the framework of three-point QCD sum rules. The cor- 
relation functions of the considered vertices when both B(D) and Kq(K)(Ki) mesons are 
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Q 2 = -m 2 B 


Q 2 = -m\ 


Average 


9B* BK 


2.40 ± 0.22 


3.62 ± 0.34 


3.01 ±0.28 



TABLE V. Value of the gs* bk coupling constant. 







A 2 


x 


Ax 


5 f D*DXi(1270)( ( 3 2 ) 


5.062 


-2.337 


1.182 


1.531 


a (D) (O 2 ) 
yDJDKltMOO^V J 


73.848 


-87.162 


118.101 


74.590 


(Xi(1270)) / 2x 
yD*DXi(1270) W ) 


0.137 


-1.507 


6.951 


1.845 


(Xi(1400)) , o2 ^ 
^0*0X1(1400)^ J 


-0.106 


1.234 


6.843 


1.847 


5 f BjBXi(1270)( ( 3 2 ) 


0.764 


0.412 


11.343 


4.708 


o (B) (O 2 ) 

y B|BXi(i4oo)vv ; 


2.463 


-2.178 


38.732 


18.980 


(Xi(1270)) , n2 N 
y BJBXi(1270)VV ) 


0.047 


-23681.595 


-31.416 


2.914 


(Xi(1400)) , n2 N 
y BJBXi(1400)VV ) 


-0.021 


0.282 


3.080 


1.233 



TABLE VI. Parameters appearing in the fit function of the coupling constants for D*DKi(1270), 
D*DKi(1400), B*BKi(1270) and B*BKi(14O0) vertices. A x and A 2 are in GeV~ l units, while 
xo and Ax are in the units of GeV 2 . 

off-shell are evaluated. The final numerical values have been obtained taking the average of 
the numerical values obtained from both off-shell cases. In the case of the axial vector Ki, 
which is either i^i(1270) or fTi(1400), the mixing between these two states have also been 
taken into account. A comparison of the obtained result on D*DK as the only previously 
calculated coupling constant among the considered strong coupling constants has also been 
made. 
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— —m D 


Cfi — -m 2 

^ ~~ ' n K 1 (mo) 


Average 


9D*DKi(1270) 


2.83 ±0.09 


1.36 ±0.14 


2.09 ±0.82 




Q 2 = -ml 


Q 2 = — 


Average 


9D*DKi(U00) 


0.97 ±0.15 


1.12 ±0.54 


1.04 ±0.78 



TABLE VII. Values of the 5'D*D_ft" 1 (i270) an d 5 , d*.da"i(1400) coupling constants in GeV 





Q 2 = -ml 


Q 2 — m 2 K ^ 1270 ^ 


Average 


5B|BXi(1270) 


1.18 ±0.07 


0.81 ± 0.45 


1.99 ±0.11 




/o2 2 

Q = -ml 


— — "^Xi(1400) 


Average 


5B s *BXi(1400) 


0.35 ± 0.05 


0.26 ± 0.04 


0.30 ± 0.05 



TABLE VIII. Values of the gB*BK!(i270) an d <?B*Bifi(i4oo) coupling constants in GeV 
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